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In a recently proposed model of "ghost condensation", spatially homogeneous states may mix, 
via tunneling, with inhomogeneous states which are somewhat similar to bubbles in the theory of 
false vacuum decay, the corresponding bubble nucleation rate being exponentially sensitive to the 
ultraviolet completion of the model. The conservation of energy and charge requires that the energy 
density is negative and the field is strongly unstable in a part of the nucleated bubble. Unlike in the 
theory of false vacuum decay, this region does not expand during subsequent real-time evolution. 
In the outer part, positive energy outgoing waves develop, which eventually form shocks. Behind 
the outgoing waves and away from the bubble center, the background settles down to its original 
value. The outcome of the entire process is thus a microscopic region of negative energy and strong 
field — "hole in the ghost condensate" — plus a collection of outgoing waves (particles of the ghost 
condensate field) carrying away finite energy. 

PACS numbers: 11.10.Lm,11.90.+t 



I. INTRODUCTION AND SUMMARY 



In view of the evidence for the accelerated expansion of the Universe, several attempts have been made recently to 
construct models in which gravity is modified at large distances PHUbJEllHISElIll O ne approach, dubbed ghost 
condensation Q , invokes a scalar field with unconventional kinetic term, like in models of k-essence @ but with the 
action depending on the derivatives only, 

S = M 4 J d 4 x P(X) (1) 

where 

x = d^d^ 

(space-time signature (+,—,—,—)). One views this model as an effective theory valid at energies below some cutoff. 
Assuming that P'(Q) < 0, excitations about a state <j> = are ghosts, hence the name of the model. A proptotype 
example of the "potential" having this property is 

P{X)= l -{X-lf (2) 

When gravity is switched off, the scalar theory is perturbatively stable provided that the background has the Lorentz- 
violating form, 

(j> = ct (3) 

where the constant c is such that 

P'(c 2 ) > (4) 
2P"{c 2 )c 2 + P'(c 2 ) > (5) 

(for the potential (J2J these inequalities imply c > 1). Cosmological evolution drives the field to a special point 
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where 

P'(cl) = 

(c* = 1 for the potential J5J). At the level of small perturbations about this background, the model has interesting 
phenomenology, as discussed in Refs. 0, El EI ■ 

There are several points to discuss beyond the perturbation theory. One is the behavior of the ghost condensate 
near sources of strong gravitational field, e.g., near black holes I n this paper we address another issue, namely, 
quantum stability of the ghost condensate. As we discuss in Section 2, even in the absence of gravitational interactions, 
there exist inhomogeneous configurations of the scalar field, bubbles, with which the state (|3J) can mix via tunneling. 
We will see that the tunneling exponent is dominated by the ultraviolet (UV) properties of the theory, so the decay 
rate is exponentially sensitive to the UV completion of the model, in accord with Ref. Q. If the UV cutoff is well 
below the scale M, the bubble nucleation rate is small. 

Still, it is of interest to understand the further evolution of the bubbles which can nucleate via tunneling. In 
conventional scalar theories with false vacua (local minima of the scalar potential), a bubble of the true vacuum, 
once created, expands practically with the speed of light Il4| , and the true vacuum eventually occupies the entire 
space. The difference between the energies of the false and true vacua is released into the bubble wall, whose energy 
thus tends to infinity at asymptotically large time. We will see in Section 3 that in the ghost condensation model 
the situation is quite different. In the central part of a bubble, the energy density is negative and, furthermore, the 
condition (0} (and also (JSJ) is violated. The system is unstable there, so the ultimate fate of the field in this region 
cannot be understood in the UV-incomplete theory. The point, however, is that this part does not expand and remains 
of microscopic size at all times; this is a "hole in the ghost condensate" . Moreover, there is no energy flow through 
the boundary of this hole 1 , so the hole has little, if any, effect on the outer part of the bubble. 

In the outer part, a shock wave (or series of shock waves) is formed 2 which propagates outwards. At the level of 
classical field theory with the action the formation of a shock wave means a singularity after which the solution 
ceases to exist. Resolving this singularity would be possible in UV complete theory only. Assuming that the singularity 
is smoothened out by the UV effects, we study the entire evolution of the bubble. We find that the ghost condensate 
behind the shock wave settles down to its original value, <p = c, = 0, so the energy of the shock wave does not 
increase in time. Therefore, the amplitude of the shock wave decreases in time; in the quantum theory the wave 
ultimately decays into quanta of the "ghost field" <f>. In the case of the special background © all these properties 
follow from the energy conservation, while for the general initial background <j) — c > we have found them by 
numerical simulations 3 . 

The outcome of the entire process is thus the hole in the ghost condensate, a microscopic region (of the size set 
by the length scale M^ 1 and/or the UV cutoff) where the energy density is negative, and the field strongly deviates 
from its background value. These holes in the ghost condensate are likely to antigravitate. 

II. NUCLEATED BUBBLES 

The energy functional of the model is 

E = J d 3 x [2P'(X)<P 2 - P(X)] (7) 

Besides the energy, there is another conserved quantity, the charge. Indeed, the field equation has the form of current 
conservation, 

d,[P'{X)d^]=Q (8) 



1 This is true, strictly speaking, only at the level of classical field theory with the action Q. Nevertheless, under an assumption that 
the energy density is bounded from below, the energy release from the hole is finite, so a possible effect of the hole on the outer part is 
minor in the quantum theory as well. 

Shock waves/kinks emerge in many models whose Lagrangians are non-linear in derivatives, see, e.g., Refs. mm. 
3 For general initial background, energy conservation does not forbid that the background behind the shock wave has lower energy density 
than the original background in front of the shock wave. Were this the case, the energy of the shock wave would increase in time, while 
the new background would eventually occupy the entire space outside the hole. The whole process would then be analogous to the false 
vacuum decay. We have not seen this phenomenon in our numerical simulations. 
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which implies the conservation of the charge 



Q 



d 3 x P'(X)<P 



(9) 



It is precisely charge conservation that renders the background stable against small perturbations: even though 
for arbitrary perturbation the energy can decrease, small perturbations in the same charge sector have larger energy 
than the background itself. 



<> (t=0, r) 



FIG. 1: An example of a configuration with vanishing energy and charge. 



<(> (t=0, r) 




FIG. 2: A bubble at the moment of nucleation. 



The latter property does not hold for large perturbations, so the state (01 may tunnel into other configurations 
which have the same energy and charge. To be specific, let us consider the background 

c = c* (10) 

(the generalization to c > c, is straightforward provided that c is not very large, otherwise one should consider trial 
configurations with non-zero spatial gradient). This configuration has 

Q = 0, E = 

so final states the system can tunnel into have zero energy and charge. An example of such a state is a spherically 
symmetric configuration with zero spatial gradient, V</> = 0, while the time derivative is a combination of step 
functions (see fig. 

4> = c\ at < r < a 
4> = C2 at a < r < b 

<j) = c* at r > b (11) 
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with < c\ < c* , C2 > . The charge and energy of this configuration are 

= a 3 P'(c 2 ) ■ ci + (b 3 - a 3 )P'(c 2 ) ■ c 2 

±-E = a 3 [2P'(c 2 )-c 2 -P(c?)] 

+ (b 3 -a 3 ){2P>(ct)-cl-P(cl)} 

The first term in Q is negative, while the second is positive (this is why the region with C2 > c* is needed). Requiring 
that Q = one finds 



b 3 -a 3 = -a 3 



P'(4) ■ c 2 



so one obtains for the energy 



where 



E = -a 3 -P'{c 2 1 )-c 1 -[F{c 2 )~F{c 1 ) 
2P'(c 2 ) • c 2 - P(c 2 ) 



F(c) 



P'{c 2 ) 



For wide class of potentials, including J2J), the function F is positive at all c, has a minimum at some point c** on the 
left of the minimum of P, i.e., < c** < c*, and diverges as c — * +0 or c — > +oo. 
As a side remark, the minimum of P occurs at the point where 

2P"(cl)cl + P'(cl) = 

It is straightforward to see that in the absence of gravity, small perturbations about the homogeneous background 
P|) grow exponentially for c** < c < c* (tachyonic region) and oscillate for c < c** and c > c* (in the former region 
they have negative energy). 

For ci < c** (that is, necessarily left of the tachyonic region) and sufficiently close to +0, there exists a point c 2 > c„ 
such that 

F(c 2 )=F( Cl ) 

With this choice of c\ and C2, both energy and charge are indeed equal to zero. The system can tunnel into such a 
configuration. Further evolution will lead to non-zero spatial gradient of 4>, because 4> is inhomogeneous in space for 
the trial configuration just constructed. Note that only the ratio a/b is determined by the energy and charge, while 
the overall spatial scale remains a free parameter; this means that the conservation of energy and charge allow for 
creation of bubbles of arbitrary size. 

It is worth noting that the step- function configuration of fig. ^ has been chosen for illustration purposes only. The 
nucleated bubbles need not have steep walls. In what follows we have in mind smooth bubble profile at the moment 
of nucleation, as shown in fig. 

The initial and final states of tunneling are somewhat unconventional here, since the field configurations are time- 
dependent. Also, when describing tunneling, one has to impose a constraint of charge conservation. The semiclassical 
formalism relevant in this situation was in fact developed some time ago [T^ . Il8j . The bottom line is that one perfoms 
the Wick rotation to the Euclidean space-time, and considers pure imaginary fields there, 

4>=-i6 (12) 

where 9 is real. The boundary conditions at initial and final Euclidean time, Tj = — oo and Tf = 0, respectively, 
require that the field be homogeneous in space, 

V6{nj) = (13) 

the latter condition being trivially satisfied for the initial state @, since at large negative times 6{x,t) = i(j)(x,t = 
—ir) = ct. Generally speaking, the semiclassical tunneling exponent equals the Euclidean action evaluated on a 
solution to the Euclidean field equations obeying the relations (|12|l and (|13|) . In the model at hand, the Euclidean 
action is still given by eq. Q), where now 

X = (d^9) 2 
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and metric is Euclidean. 

It is straightforward to see, however, that the model as it stands does not have relevant Euclidean solutions. 
Indeed, under rescaling 

6(x>*) -> A0(A" V) 

the boundary conditions (|13f) and the initial configuration 6 = ct do not change, while the action scales as 

A 4 S 

Thus, there are no relevant saddle points; configurations of small size have small Euclidean action. 

The scaling argument tells that the bubble nucleation rate is UV dominated. If the UV cutoff A is smaller than 
M, the decay exponent is large, 

M x 4 



S~[ T ) (14) 



so the decay rate is small 4 , in accord with the analysis in Ref. Q. 



III. REAL TIME EVOLUTION 
A. Classical evolution in (l+l)-dimensions 

1. Legendre transformation 

To understand how configurations of the type shown in fig. [21 evolve in Minkowski space-time, let us first neglect 
the spatial curvature of the bubble walls, and thus consider the system with the action (Q in (f + I) dimensions. The 
explicit form of the field equation (JSJ in (1 + I) dimensions is 

(P' + 2P"4> 2 t )cp, tt - AP"^ x cf> >t( f> !Xt - (P' - 2P"4,%)(t>, xx = (15) 

Here and in the following comma denotes differentiation. This non-linear equation is simplified by the Legendre 
transformation. Instead of x, t, one chooses new independent variables 

Z = <t>,t 

V = 0.x (16) 

and, instead of <j), one considers an unknown function 

v) = -<t> + + x v 

This change of variables is legitimate in any region where the pairs of coordinates (£, rf) and (x, t) are in one-to-one 
correspondence. Equivalently, in such a region the Jacobian of the coordinate transformation, 

J = U££U tVV - (u^) 2 (17) 

is nowhere zero or infinite. The second derivatives are 



(18) 



>,« 


= J- 1 




\xt 


= -J- 


1 ■ u Xv 


,xx 


= J- 1 





4 One way to understand this property is to consider higher order operators added to the action Q. With a single additional operator of 

4ti — 12 

the form A - ( 2n-6 ' (d™<f>) 2 with n > 3, the same scaling argument gives the estimate S ~ (M/A) for the action at a saddle point. 

Adding all higher order operators effectively corresponds to the limit n — > oo, in which the estimate 1141 is obtained. 
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while the original coordinates arc related to £ and rj as follows, 

t = u x (19) 
x = u n (20) 

The original equation (|15J) in terms of new variables takes the form 

(P' + 2P"£ 2 K W + 4P"^„ - (P' - 2P'V K« = (21) 

where P is a function of the combination (£ 2 — rj 2 ). This is now a linear equation with coefficients that depend on 
the coordinates £, 77. Once a solution u(^,rj) is found, the quantities of interest <fij = £(2^) and <p tX = rj(x,t) are 
obtained, at least in principle, by inverting eqs. I|19fl . (|20|l . 

To simplify eq. I|21(l , let us take advantage of its Lorentz invariance and introduce new variables p and 9 related to 
£ and 77 as follows, 

t = 4>, t = P cosh 9 (22) 
ri = <t>, x = psinh(9 (23) 

(of course, we made an assumption here that \(f>^\ > \<l>,x\'i it is valid in all cases we consider). Clearly, the meaning 
of p is that it is a scalar 



p = J{d^Y = 



In terms of the variables p and 9, the original coordinates are 

1 



Now eq. l|2~T|) reads 



t = u „ cosh 9 u sinh 9 (24) 

P ' 

x = — u,pSinh#H — ugcoshfl (25) 



p' + spy p , p + 2py 

— 2 u ,ee - P u^ P p Up = (26) 



where P = P(p 2 ), and prime still denotes the derivative with respect to its argument p 2 . 

This equation is hyperbolic when P' and (P' + 2P"p 2 ) have the same sign, i.e., when the small perturbations of 
the field (f> do not grow (no tachyons). This equation is elliptic in the tachyonic region, and has singularities (zeroes 
in front of one of the second derivative terms) on the boundary of that region, i.e., at P' = and at P' + IP" p 2 = 0. 
The second boundary point is not of interest for our purposes; about the first boundary point we will have to say 
more later. 

Let us consider the region on the right of the minimum of P, where both P' and P" are positive. The final 
simplification of eq. (|26|) is made by introducing a new variable r instead of p, 

(27) 




pi pi 

where for future convenience the constant of integration is chosen in such a way that 

r(p = c») = 



Then instead of eq. H26J) one obtains 



where 



u.ee — 7i, Tr — F{t)u, t = (28) 
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It is instructive to point out that near the point c* one has 

P' = const • (p 2 - cl) 
Then the above expressions have simple form. In particular, 



and eq. reads 



r = 2J P —^ (29 ) 



u 00 — u TT u T = (30) 

r ' 



The latter equation has the form of the wave equation in (2 + 1) dimensions for 0(2)-symmetric functions. 



2. Shock waves 



Let us consider initial data such that 4>,x = everywhere (cf. eq. \i'6\ ) and <fij > c*, i.e., P' and P" are positive 
on an entire line — oo < x < +oo. Let us further assume that at t = is a monotonically increasing function of x, 
which tends to certain finite values as x — > ±oo, 

<j) t t{t = 0, x — > — oo ) — > p- (31) 
<t>, t {t = 0,x -> +oo) -> p + (32) 

with 

P+> P-> c* 

This initial configuration is shown in fig. El and it corrsponds to the outer part of the bubble, namely, the region r > b 
for the configuration of fig. [21 To this end, x is to be identified with — r (thus, the bubble is meant to be on the right 
in fig. |21 this somewhat bizarre convention simplfics the further discussion) . 




FIG. 3: Initial data in the outer region of a bubble. 



These initial data are translated into the Legendre variables as follows. According to eqs. (|23(l . (|24(l and l|25|l . the 
initial condition <fi. x = is satisfied when the line t — coincides with the line 9 = 0, and u = const on that line. 
Thus, the initial data are 

u, e + , 6 = (33) 
Up = , 6 = (34) 

Once the initial data <fi,t(t = 0,x) is specified, the initial data for the Legendre variables is obtained by inverting the 
relations 

p = <f> t t(t = 0,x) 

x = -u e (6 = 0,p) (35) 
P ' 
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According to eqs. (12211 and l|31(l . i|32[l . the initial data for u,g are specified on a finite interval 

9- < P < P+ (36) 

Since <j> t at t = is an increasing function of x, eq. H35|) means that u t g at = is an increasing function of p. Since 
x runs from — oo to +oo, cq. i|35|) implies that u t g is singular at the ends of the interval (|36[) . 

Ufi{6 = 0, p — > p-) — > — oo (37) 
u, 9 (0 = 0,p^ p+) -> +oo (38) 

Instead of p, one can of course use the coordinate r in the above expressions. 






Let us now discuss the behavior of the solution u(t,8) to cq. 128|) . in the (t, 0) plane. The singularities (|37|l and 
(13811 that start at r_ and t + , respectively, move along the light-like lines on this plane. The relevant light-like lines 
are the lines L and R in fig. 01 (the other two lines are screened by the lines L and R). 

The right singularity, labeled R in fig. ^ occurs at 

6 = t + -t (39) 

Near the right singularity the function u has the form 

u = A{t)C+(6 + t) (40) 
where A(t) is a slowly varying function, while C' + (z = (6 + r)) varies rapidly and has the property 

C+ (z) -> +oo , z -> r+ - 

The form of C+ is determined by the singular part of the initial configuration u_g{6 = 0, p) as p — > p + , i.e., by the 
behavior of </>,t(t = 0, x) as x — > +oo. 

^From eqs. I|24|l and l)25|l one finds that near the right singularity 



W) 

Mr) 
P(r) 



t = -}~lc' + (e + t)( P t p cosh9 - sinhO) (41) 
C' + (6 + T)(-pT p smh8 + cosh8) (42) 



Note that from eq. I|27|l it follows that pr p > I, so time t is positive near the right singularity. The coordinate x is 
also positive near the right singularity, at least in the lower right corner of fig. where 9 is small. Thus, the lower 
right part of the line R corresponds to the asymptotics t — > oo, x — > +oo. 



9 




FIG. 5: Line of constant (<9 M 0) 2 in (r, 6>)-plane. 



A line of constant r, i.e., a line of constant X = (dp4>) 2 , shown in fig. [5] hits the singularity. According to eqs. 141|) 
and H42JI . at large t and x this line becomes a straight line in (x, i)-plane, with the slope 

cos1i(t + — t) — pr p sinh(r + — r) 



x 
1 



pr^p cosh(r + — r) — sinh(r + — r) 



(43) 



where we made use of the relation <|39[) to express through r. For r sufficiently close to r+, the lines of constant X 
move towards positive x, so there is a wave moving right. In the case of the bubble of fig. this wave moves inwards 
the bubble. 

For generic initial data there is a point on the right singularity line R where the two terms in the numerator in 
eq. I)43|) cancel each other (shown by a circle in fig. 0. Then left of this point (i.e., at t < to) the coordinate x 
is negative (and tends to — oo as the line of constant r approaches the singularity line). In this region, the lines of 
constant (dptfi) 2 move left, so there is a left-moving wave. In between the left-moving and right-moving waves, the 
space is gradually occupied by "constant" ghost condensate 



i>.t = p(ro) cosh(r + - r ) 
b,x = p(t~o) sinh(r + - r ) 



(44) 



The final point to mention about the right singularity line is that the slope (|43|l increases as p increases, so the lines 
of constant X do not intersect on (x, t)-plane, sec fig. Indeed, one finds 



2p 3 



^l + 2P"p 2 /P' 



[pT tP cosh(r+ - 
pin pin 



t) — sinh(r + — r)] 5 



P 
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P> 



(45) 



which is positive for large class of the potentials P. As it evolves, the wave gets spread over larger interval of x. The 
asymptotic velocity of the wave is always smaller than the speed of light: the maximum slope (|43f) is attained in the 
lower right corner where r w t + and 8 w 0, and there 



x 
1 



1 



O t ,p)0+) 



< 1 



For initial data close to the minimum of P, the maximum slope, and hence the velocity of the wave arc small, since 
pr p is large, see eq. I|29|). 

Let us turn to the left singularity line L. This line corresponds to 



6 = t - T_ 
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and the solution near the singularity line has the form. cf. cq. (|40|) . 

u = A(t)C-(t-0) (46) 
The function C_ is rapidly varying; it is an increasing function of its argument with the property 

C_ (z) -> -oo , z = (t - 0) -> r_ + 

At the very first sight the situation here is similar to that near the line R. This is not quite true, however. The 
expressions for t and x near the line L (analogs of eqs. 1)41(1 and 1)42(1*1 arc 

i = 4t c -( t - 8 )(p t ,p c ° sh 61 + smh °) ( 47 ) 



r 



z = -^hlc'_(T-9)(pT p smh8 + cosh6) (48) 

Time t is again positive while x is now negative. Thus the left part of fig. 0] describes a wave moving left. In the 
context of a bubble of fig- El this wave moves outwards. 

The point, however, is that the absolute values of the slopes of lines of constant r (i.e., constant X) on the (x,t)- 
plane increase as p increases: the relevant formula is again given by eq. 1(45(1 . but with the opposite overall sign (the 
wave moves left) and with (r — t_) substituted for (r+ — r). Thus, at given large t, points (in x-space) with larger p 
are to the left of points with smaller p. On the other hand, at t = the situation is opposite, since the initial data are 
such that (f> tt — p increases with x as shown in fig. |3J This means that lines of constant X intersect on (x, t)-planc, 
see fig. Of course, the whole treatment breaks down when these lines intersect, as there emerges a singularity on 
(x, t)-plane. This singularity is of the type of shock wave, or kink, at which the first derivatives and <f> tX are step 
functions: just before the lines of constant X intersect, the values of X(x), and hence 4> t t and 4>,x are substantially 
different at neighboring points of space. From the above analysis it follows that this shock wave moves left in fig. |3 
at least at the moment when it gets formed — in the context of a bubble this is the motion outwards. 



3. Numerical analysis 



Equation 115|) for the evolution of the ghost field can be studied numerically. For the purpose, we found it convenient 
to take £ = <j> t and 77 = (f) iX as new dependent field variables 5 . In terms of these, the equation of motion can be recast 



Here we consider £ and r) still as dependent functions of t and x, i.e. we do not perform the Legendre transformation of subsection 

imm 
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into the form 



(P' + 2P"e)^t - 4P"r)ti,x - (P' - 2P'V k* = 



where P = P(£, 2 — i] 2 ), with the additional constraint 



lit = L 



(49) 



(50) 



We discretize now the space dependence of the variables £, 77 by introducing a lattice with uniform lattice spacing a 
and defining 



r) n = 77(710) 



(51) 



We also make the overall lattice finite imposing — N < n < N , introduce Neumann boundary conditions on </>, and 
replace the derivatives ^, x ,?7,x with their central difference approximations (V£) n = (£ n +i — £n-i)/(2o), (Vv)n = 
{Vn+i — ?7n-i)/(2a). After these steps eqs. Q49fl. I|50|) take the form of a set of coupled ordinary differential equations 
of the first order for the time dependence of the variables £, n ,Vn- 



rtt 
drj 
It 



4P"t7^VC + (P' - 2P'V )Vt7 



(p/ + 2P"C 2 ) 



(52) 
(53) 



We integrated these equations in time by using the second order Runge-Kutta formula. 

We would like to make a few observations about our numerical procedure. The main concern in integrating evolution 
equations such as the ones at hand is the possible onset of numerical instabilities. The avoidance of such instabilities 
limits the maximum time step to values that are generally small enough to render the use of higher order integration 
formulae unwarranted. This is why we used the second order Runge-Kutta algorithm. The adequacy of this algorithm 
was confirmed by the conservation of energy, which, in the stability region, was accurate to better that one part per 
million. However, the time integration can become unstable. This happens a) when the formation of a shock wave 
generates a wave-front spanning only a few lattice spacings and b) when the signature of the equation is elliptic over 
some domain of values of x. In the case (a) the approximation to the space derivative obviously fails, and this leads to 
instabilities in the time integration. We have noticed, incidentally, that the onset of such instabilities is more serious if 
one uses <fi itself, rather than 77 = <p >x , as one of the independent variables. The use of rj avoids the need of discretizing 
second derivatives with respect to x and leads to a better behaved integration algorithm. From the physical point of 
view, one can also argue that (f> tX plays a more fundamental role than <fi in the system under consideration. 

The instabilities that occur in the case (b), namely in the domains where the equation is elliptic, are of physical 
nature, and not caused by shortcomings of the integration procedure. Indeed, with an elliptic equation, perturbations 
grow exponentially when one integrates the equation forward in time. The damping of such instabilities would be left 
to the UV completion of the theory. 

Whether due to a shock wave or to a domain of elliptic signature, numerical instabilities are characterized by an 
exponential growth of the short range fluctuations. In order to tame such instabilities we supplemented our evolution 
algorithm with an additional step, which serves to dampen the growth of those fluctuations. Namely, after the 
standard Runge-Kutta integration step, we use a fast Fourier transform (FFT) to calculate 



Vk 



(54) 



and dampen £, f\ by replacing 



Vk 



£/ _ e -a[l-cos(27rfc/iV)] dt/a 2 ^ 
fj> = e -a[l-co B (2irfe/JV)] dt/a 2 ~ fe 



(55) 



dt being the integration time step and a suitable parameter. After this step, we go back to the original variables by 
the inverse FFT. 
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We have found that, with a judicious choice of the damping factor a, we were able to evolve the equations well 
past the formation of the shock wave and also for some time in the presence of an elliptic domain, while preserving 
the conservation of energy to a few per cent, or better. Of course, the front of the shock gets slightly smoothed out 
and, in the elliptic domain, the short range fluctuations, after some initial wild growth, get gradually damped down 
and become intolerable only at relatively late times. It is interesting to observe that our procedure, while mainly 
motivated by the desire to keep the numerical integration under control, could also be thought of as a way to mimic 
the effects of the UV completion. 

Finally, the discretization of the equations for the three-dimensional evolution with spherical symmetry proceeds 
along the same lines, after having inserted the appropriate metric factors (powers of r). 

In fig. we show the profile of the field configuration obtained by solving eq. (|15fl numerically, with the potential P 
given by eq. J2J) and initial data similar to those shown in fig. [21 This solution clearly demonstrates the formation of 
the shock wave and its motion to the left. There is also a soft wave slowly propagating to the right. In our numerical 
simulations we always observed this behavior of solutions for initial data of the type shown in fig. [SJ 
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4- Around the critical point 

Let us now consider another initial configuration, again with <f> x = but with the profile of cf> t shown in fig. [S] 
This configuration corresponds to the region r < a in fig. [21 i.e., the region near the point where 4>,t — c * ( m fig- [HI 
the bubble is meant to be on the left, like in fig. El and unlike in fig. |3J. At the left of the critical point t = c* the 
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system is unstable. Let us study whether this instability propagates into the region on the right of this critical point, 
and also study the motion of this critical point as time evolves. 




FIG. 8: Initial data in the interior of a bubble. 



Let us consider the behavior of the solution <j){x, t) in the region where X > cj. In this region one can still use the 
variables p and 9 (or, equivalently, r and 9) in the Legendre conjugate problem. Now, the critical point X = p 2 = c 2 
corresponds to r = 0, and in the vicinity of this point the function u(r, (9) obeys eq. (|3(J[) . ^From the analogy to 
0(2)-symmetric (2 + l)-dimensional problem it is clear that the solution and its derivatives never become singular 
at t = 0, provided initial data are smooth. Futhcrmore, the solution is uniquely determined by the initial data at 
X > c 2 .. This means that the behavior of the solution <p{x,t) in the region X > c 2 is not sensitive to its behavour 
at X < c 2 ; the regions right of the critical point and left of the critical point do not talk to each other; whatever 
happens in the inner, tachyonic region of the bubble, has no effect on the outer region 6 . 

Let us see that the critical point X = cl moves left in fig. [SJ i.e., inwards the bubble. To this end, we note that 
the initial data for the Legendre problem arc again formulated in a finite interval of p, which is now c* < p < p+, 
i.e., in a finite interval < r < r+. The initial data is u(9 = 0, r) = 0, while u^(9 = 0, r) is singular at r = r + . At 
non-zero 9, the singularity propagates along the light-like line 9 = t+ — r, as shown in fig. El Near the singularity line 
the solution still has the form given by eq. 1)40(1 . and this line again corresponds to t — > oo. The slope ((43JI is positive 
on the lower right part of this line where (r + — r) is small, so the wave corresponding to the region near the corner 
9 = 0, t = t + is again moving right in fig. |SJ This is now the motion outwards the bubble. 



e 




T=0 X 



FIG. 9: Singularity line in (r, #)-plane for initial data of fig.|Hl 



In the classical field theory problem, a singularity may eventually be formed in the tachyonic region, and the solution does not globally 
exist at later times. 
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At p close to c* one has pr p = {p/c* — l) -1 / 2 , sec ca. <|29[) . which is large. Thus there is always a point denoted 
by a circle on the singularity line R, at which the slope (|43|) vanishes. Regions right and left of this point correspond 
to motion right and left in fig. |H1 i.e. outwards and inwards the bubble, respectively In particular, the critical point 
r = (i.e., X = c 2 ) moves left in fig. 00 i.e., inwards the bubble. 

Another, direct way to see that the regions with X < c 2 and X > 
consider eq. Ijl5(l itself. Rescaling the field in such a way that 



1 



one has near the critical point 



P' = X - 1 

Equation Ijl5|) then reads 

(3$ - <t>% - l)4>,tt ~ H,^,t4>,xt - {4> 2 t ~ 34>% - = (56) 

This may be viewed as a wave equation with coefficients depending on <j)^ x and <f>^. Let us find the characteristics of 
this equation that starts at a point (x,t) where X = <f> 2 t — <p x = 1, i.e., at the critical point. At this point, eq. I|56|) 
reduces to 

§ 2 t4>,tt - 2^,x0,t0,^t + 4> 2 x 4>,xx = (57) 
Thus, the equation for the characteristics is 

(j> 2 x ■ dt 2 + 2<j> tX </> )t ■ dxdt + 4> 2 t ■ dx 2 = 
We see that the characteristics are degenerate and obey 

dx d> T / \ 

Tt=-t m 

Let us now see that the critical point where X = 1 moves precisely along this characteristic. The motion of the critical 
point is determined by the equation 

d(0 2 t - <!>% - 1) = ($ - <f>% - l), t dt + - 4> 2 X - l), x dx = 

which gives 



'>,x<P,xt 



dx _ <t>,t<t>,t 
dt 4>,t4>,xt - 4>,x4>,xx 

Making use of eq. (|57|l we find that the right hand side of this expression is precisely the same as the right hand side 
of eq. (JSHJl, so the charactristic and the world line of the point where X = 1, indeed coincide. 

Since the line X = 1 is a characteristic, signals emitted on the left of this line do not propagate to the right of this 
line, and vice versa, i.e., the regions with X < c 2 and X > c 2 indeed evolve independently. Equation (|58(l determines 
the motion of the critical point; for an initial configuration shown in fig. [S] the velocity of this point initially vanishes 
(since 4>^ x = at i = 0), and then becomes negative, as positive 4>,x develops. The critical point indeed moves inwards 
the bubble. 

One way to understand the unusual property that the evolution in the regions with X < c 2 and X > c 2 proceeds 
independently, is to recall the expressions for the energy- momentum tensor and conserved current, 

= 2P'{X)d ll <j>d v <j>-P{X)r lliV (59) 
j„ = P'(X)d^ (60) 

One observes that all components of T^ v and are equal to zero at the surface in space where X = c 2 . Thus, there 
is no transfer of energy and charge between the regions with X < c 2 and X > c 2 . This explains, at least partially, 
the independence of the evolution in these two regions. It is worth stressing that this argument works in any number 
of dimensions. 
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FIG. 10: Initial stages of evolution of a bubble in (1 + 1) dimensions. Profiles of X are shown at different times. The irregular 
pattern at x < —40 for t — 22 is due to instabilities in the tachyonic region where X < c\ = 1. 



5. Overall evolution 

Combining the two pictures, corresponding to the initial data of figs. O and |H1 we obtain the following qualitative 
properties of the evolution of a bubble, still in (1 + 1) dimensions. The inner part of the bubble is unstable due to the 
tachyonic character of the region < X < c%. Nevertheless, the region X > is not sensitive to this instability. 
The critical point where X = moves towards the center of the bubble, so the region of large fields (hole in the ghost 
condensate) remains of small size. In the outer part, a wave moving outwards the bubble develops. Eventually this 
wave forms a kink. We show in fig. 1101 the evolution of a bubble, which we obtained by solving eq. I|15fl numerically 
with the bubble- like initial data (in particular, initially <f) tX = 0). The region left of x ~ —20 in this hgure (units of x 
and t are arbitrary) has initially X < 1. The field at large negative x is thus unstable; this instability is manifest in 
fig. 1101 The motion of the boundary of the tachyonic region to the left (inwards the bubble) and a shock wave moving 
right are also clearly visible. 

B. Four dimensions 

Although the above analysis has been performed in (1 + 1) dimensions, we argue that our main findings are valid 
in (3 + l)-dimensional theory as well, at least for 0(3)-symmetric bubbles. At large distances from the bubble center, 
the curvature of a sphere r = const is small, so the evolution is similar to that in (1 + l)-dimcnsional theory. There 
are two consequences of this simple observation. First, the tachyonic region, where X < ej, does not expand; instead, 
it shrinks and forms a hole in the ghost condensate. The energy of this region is negative but finite at the moment of 
nucleation. According to eq. (|59|l . there is no transfer of energy from the hole to the outer region where X > cjjj, so 
the overall energy of the outer region, referenced from the energy of the original background <j> = c, remains a finite 
constant 7 . 

Second, an outgoing wave of positive energy eventually forms a shock. It is conceivable that the UV effects smoothen 
out the profile of the wave. In front of the wave, i.e., at large r, the ghost field is still in its original state, V</> = 0, 
4> = c. Let us discuss the field behind the wave. 



7 Assuming that the UV effects render the energy of a microscopic hole finite, the latter property is independent of the form H59I of the 
energy-momentum tensor, which is not exact in a complete theory. 
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Let us first consider the special value of the original background, 4> = c*. In that case the conservation of energy 
requires that the background behind the outgoing wave is also <f> = c*, V0 = 0. This follows from the fact that the 
energy density of the original background is zero, while the overall energy of the positive energy region (outside the 
hole) is finite. This implies that X = c 2 behind the outgoing wave, otherwise the energy of the positive energy region 
would increase as r 3 (t) where r(t) is the radius of the wave, which grows in time, r(t) — > oo as t — > oo. The property 
that X = c 2 behind the wave, does not guarantee by itself that the background behind the wave is the same as the 
original background. However, the energy density of the wave decreases in time (otherwise its total energy would 
grow as r 2 (t)). Therefore the amplitude of the wave decays, which is only possible if V</> = 0, cb = c* both in front of 
and behind the wave 8 . 

The above argument does not work for the general values of the original background, <j> = c > c*. In that case 
the conservation of energy and charge does not forbid that the energy density behind the outgoing wave be smaller 
than the energy density of the original background (i.e., of the field in front of the wave). However, in our numerical 
simulations the background behind the outgoing wave(s) always settled down to its original value, so the energy of 
the wave did not grow in time, and the amplitude of the wave decreased as it moved away from the bubble center. 
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FIG. 11: Initial stages of evolution of a bubble in (3 + 1) dimensions. Profiles of y/~X are shown at different times. Outgoing 
and incoming waves are formed, the incoming wave eventually bounces from the origin. 

In numerical simulations, we observed a fairly complex behavior of the system even in the 0(3)-symmetric case. As 
a reservation, we did not incorporate the region where X < c 2 into our numerical study of the long-time properties 
of the solutions. The reason is that the system is unstable for X < c%, and this instability becomes intolerable in 
numerical simulations too early. Instead, we considered initial data such that Vc5 = 0, > c*. To mimic possible 
effects of the boundary of the stable region (a surface at which X = cl ) , we considered various boundary conditions 
near the origin, namely (i) smoothness at the origin, <p_ r = at r = 0; (ii) free boundary condition at some fixed 
r = rh (meant to be the radius of the hole); (iii) X = (j> 2 t — <fi 2 r = c 2 at r = r^. The overall behavior of the solutions 
at r r/j was essentially independent of the choice of the boundary condition. 

A typical sphericaly symmetric solution in the theory with the potential (0) is shown in figs. II II 1141 The initial data 
for this solution are such that <^ r = while 4> has the bubblcdike shape, as shown in fig. ^2 The initial stages of the 
evolution, shown in figs. II II — 1131 are indeed quite complex, but at late times the solution is simple: there is a sequence 
of outgoing waves with the backgrond behind them equal to the background in front, see fig. 1141 The amplitudes 



8 This argument does not work in (1 + 1) dimensions: the conservation of energy and charge do not forbid that the background behind the 
outgoing wave is Lorentz-boosted with respect to the initial background. We have not observed such a phenomenon in our numerical 
simulations. 
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FIG. 12: The second outgoing wave forms a kink. Note that the value of X near the origin temporarily becomes smaller than 
in the original background. 
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FIG. 13: The second kink moves back and hits the origin. 



of the outgoing waves indeed decrease in time, so the system settles down back to the homogeneous configuration 

4>,r = 0, </> = C. 

Thus, the outcome of the entire process discussed in this paper is a hole in ghost condensate plus a few outgoing 
"ghost waves" (in quantum theory the latter are particles of the ghost field <j>). Without knowing the UV-complcte 
theory one cannot tell much about the properties of the negative energy holes. Our analysis implies that these holes 
do not expand, and thus remain of microscopic size. Although in the theory with the action Q per se the field inside 
the hole is unstable, it is conceivable that the state of the hole gets stabilized due to UV effects. It is of interest to 
further investigate the phenomenology of ghost condensate with holes. 

The authors are indebted to N. Arkani-Hamed, F. Bezrukov, S. Dubovsky, D. Gorbunov, A. Gorsky, D. Levkov, 
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FIG. 14: The third outgoing wave is formed; the background behind the outgoing waves settles down to its original value. 
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